Polymer heat transport enhancement in thermal convection: the case of 

Rayleigh- Taylor turbulence 
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We study the effects of polymer additives on turbulence generated by the ubiquitous Rayleigh- 
Taylor instability. Numerical simulations of complete viscoelastic models provide clear evidence that 
the heat transport is enhanced up to 50% with respect to the Newtonian case. This phenomenon is 
accompanied by a speed up of the mixing layer growth. We give a phenomenological interpretation 
of these results based on small-scale turbulent reduction induced by polymers. 

PACS numbers: 



Controlling transport properties in a turbulent flow is 
an issue of paramount importance in a variety of situa- 
tions ranging from pure sciences to technological applica- 
tions [l|-|3|. After Toms [4|], one of the most spectacular 
way to achieve this goal consists in adding inside the fluid 
solvent a small amount of long-chain polymers (parts per 
million by weight). The resulting fluid solution acquires 
a non-Newtonian character and the most interesting dy- 
namical effect played by polymers is encoded in the drag 
coefficient, a dimensionless measure of the power needed 
to maintain a given throughput in a pipe. With respect 
the Newtonian case (i.e., in the absence of polymers), it 
can be reduced up to 80% H@|. 

In many relevant situations {e.g. atmospheric convec- 
tion) the velocity field is two-way coupled to the tem- 
perature field with the result that, together with mass, 
also heat is transported by the flow. Because drag re- 
duction is associated to mass transport enhancement, an 
intriguing question is on whether this is accompanied by 
a similar variation in the heat transport. 

In this Letter we demonstrate the simultaneous occur- 
rence of mass transport enhancement (drag reduction) 
and heat transport enhancement induced by polymers 
in a three-dimensional buoyancy driven turbulent flow 
originated by the ubiquitous Rayleigh-Taylor (RT) in- 
stability 0, HJ. This instability arises at the interface 
between a layer of light fluid and a layer of heavy fluid 
placed above and develops in a turbulent mixing layer 
(see Fig. [TJ which grows accelerated in time. Heuris- 
tically, the RT system can be assimilated to a channel 
inside which vertical motion of thermal plumes is main- 
tained by the available potential energy. Our idea on 
the possibility of observing drag reduction in this system 
is suggested by recent analytical results which show a 
speed-up of the instability due to polymer additives 
Moreover, examples of turbulent drag reduction without 
boundaries have been recently provided, e.g., in [Iol - [l3| 

Direct numerical simulations of primitive equations 
show that thermal plumes are faster in the presence of 
polymers (see Fig. [TJ, therefore the mixing layer acceler- 



ates (up to 30% at final observation time) with respect 
to the Newtonian case and complete mixing is achieved 
in a shorter time. A second and more dramatic effect, 
also clearly detectable in Fig. [TJ is that polymers reduce 
small scale turbulence |10h12| . As a consequence, ther- 
mal plumes in the viscoelastic case are more coherent 
and transport heat more efficiently. Quantitatively, the 
enhancement of the heat transport corresponds to larger 
values (more than 50% at final observation time) of the 
Nusselt number with respect the Newtonian case. 




FIG. 1: Vertical sections of temperature field for Newto- 
nian (left) and viscoelastic (right) RT simulation at time 
t = 2r starting from the same initial conditions. White 
(black) regions correspond to hot (cold) fluid. Boussinesq- 
Oldroyd-B equations (TJJ are integrated by a standard fully 
dealiased pseudo-spectral code on uniform grid at resolution 
512 x 512 x 1024. Physical parameters are Pr — vjn — 1, 
n — 0.2 (?7 = for Newtonian run), fig — 0.5, 8q = 1 
(Ag = 0.25). Deborah number De = t p /t is De = 0.2. The 
initial perturbation is seeded in both cases by adding a 10% 
of white noise (same realization for both runs) to the initial 
temperature profile in a small layer around the middle plane 
= 0. 

We consider the incompressible RT system in the 
Boussinesq approximation generalized to a viscoelastic 
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fluid using the standard Oldroyd-B model 14 j 
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together with the incompressibility condition V ■ u = 0. 
In (JlJ T(x, t) is the temperature field, proportional to 
the density via the thermal expansion coefficient (3 as 
p = po[l — f3(T — To)] (po and To are reference values), 
Gij (x, t) is the positive symmetric conformation tensor of 
polymer molecules, g = (0, 0, — g) is gravity acceleration, 
v is the kinematic viscosity, n is the thermal diffusiv- 
ity, rj is the zero-shear polymer contribution to viscosity 
(proportional to polymer concentration) and r p is the 
(longest) polymer relaxation time (l4j |. 

The initial condition for the RT problem is an unstable 
temperature jump T(x, 0) = — (9o/2)sgn(z) in a fluid at 
rest u(x, 0) = and coiled polymers <r(x, 0) = I. The 
physical assumptions under which the set of equations 
([1]) is valid are of small Atwood number A = (l/2)/3(9o 
(dimensionless density fluctuations) and dilute polymer 
solution. Experimentally, density fluctuations can also be 
obtained by some additives (e.g., salts) instead of temper- 
ature fluctuations: within the validity of Boussinesq ap- 
proximation, these situations are described by the same 
set of equations ([1]). In the following, all physical quanti- 
ties are made dimensionless using the vertical side, L z , of 
the computational domain, the temperature jump 9q and 
the characteristic time r = {L z / Ag) 1 / 2 as fundamental 
units. Elasticity of the polymer solution is measured by 
the Deborah number De, the ratio of polymer relaxation 
time to a characteristic time of the flow. In our unsteady 
case De grows in time starting from De ~ 0, therefore 
viscoelastic effects are initially absent. An estimate of 
the largest Deborah number achievable is based on the 
large scale convective time as De = t p /t. 

Total energy of the solution has an additional elastic 
contribution to kinetic energy E = K + £ = (1/2) (it 2 ) + 
(ur}/T p )(tra) and the energy balance for ^ reads 

dP „ , „, dE 2vn .. 

— = Pg{wT) = — +e„ + -^ [{tra) - 3] (2) 



dt 



dt 



where P — —f3g{zT) is the potential energy and £ v = 
v((d a up) 2 ) is the viscous dissipation and the last term 
represents elastic dissipation. Because this last term in 
(|2|) is not negative, one might expect that the presence of 
polymers accelerates the consumption of potential energy 
with respect to the Newtonian case (77 = 0), as it is indeed 
observed in Fig. [5J 

Of course, the speed-up of potential energy consump- 
tion due to polymers does not automatically imply the in- 
crease of kinetic energy growth. Part of potential energy 
is indeed converted to elastic energy £ by polymers elon- 
gation. The inset of Fig. [5] shows indeed that kinetic en- 
ergy for viscoelastic runs is larger than in the Newtonian 
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FIG. 2: Time evolution of total energy E = K + S, kinetic 
energy K (upper inset) and potential energy P (lower inset) 
for the Newtonian run (De = 0, black continuous line) and 
the viscoelastic run {De = 0.2, red dashed line). 



case (of about 40% at t — 2.5r). This is the fingerprint of 
a "drag reduction" as defined for homogeneous-isotropic 
turbulence in the absence of a mean flow 0, EH. 
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FIG. 3: Time evolution of velocity anisotropy w r , ns /u rms 
where w and u are the vertical and one horizontal veloc- 
ity components, respectively. Black line is the Newtonian 
run at De — 0, red dashed line is the viscoelastic run at 
De — 0.2. Inset: the evolution of the ratio of velocity gradi- 
ents (dzw) rrns /(d x u) rms . 

The most important effect of polymers on turbulent 
velocity is to generate more coherent thermal plume with 
respect the Newtonian case, as it is evident in Fig.[T] This 
reflects in larger vertical component of the velocity with 
respect the horizontal one, i.e., an increased anisotropy of 
the velocity field. This effect is evident in Fig. [3] where we 
plot the ratio of vertical rms velocity w rm s to horizontal 
one u r ms- The anisotropy ratio, which is around 1.8 for 
the Newtonian case [l5|], becomes larger than 2.5 for the 
viscoelastic run. More important, in the viscoelastic case 
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the anisotropy persists also at small scales (i.e., in the 
ratio of velocity gradients), while it is almost absent in 
the Newtonian case (see inset of Fig. [3]). 




FIG. 4: Growth of the mixing layer thickness ho.gs(t) de- 
fined as the vertical range for which |T"(«)| < 0.98#o/2 as a 
function of dimensionless time t/r for Newtonian run (black 
line) and the viscoelastic run (dashed red line) starting from 
the same initial condition. Upper inset: mean temperature 
profiles T(z) for the two cases at time t = It. Lower in- 
set: temperature variance profiles a%(£) at different times vs 
z/h(t) (black Newtonian, red: viscoelastic run). 

Despite the fact that RT turbulence has vanishing 
mean flow, a natural mean velocity is provided by the 
growth of the width h(t) of the turbulent mixing layer 
where heavy and light fluids are well mixed. For ordi- 
nary fluids at small viscosity, as a consequence of con- 
stant acceleration, one expects h(t) — aAgt 2 where a is 
a dimensionless parameter to be determined empirically 
l3-18|- Several definitions of h(t) have been proposed, 
based on either local or global properties of the mean 
temperature profile T(z,t) (the overbar indicates aver- 
age over the horizontal directions) fljl- 22|. The simplest 
measure h r is based on the threshold value of z at which 
T(z,t) reaches a fraction r of the maximum value i.e. 
T(±h r (t)/2,t) = Tr6 /2. 

Figure 0] shows the growth of the mixing layer thick- 
ness for both Newtonian and viscoelastic RT turbulence. 
As already suggested by Fig. [TJ in the viscoelastic solu- 
tion the growth of the mixing layer is faster than in the 
Newtonian case (i.e. larger h(t), up to 30% at t = 2.5r), 
therefore we have an effect of polymer drag reduction, 
i.e. polymer addiction makes the transfer of mass more 
efficient. The inset of Fig. 2] shows that the increased 
efficiency is a global property of the mixing layer and 
the temperature profile of the viscoelastic solution cor- 
responds to the profile of the Newtonian case at a later 
time. 

Also in Fig. |4] we plot the variance profiles of temper- 
ature field computed at different times for both Newto- 
nian and viscoelastic turbulence. In both cases, o~\(z) 



at different times collapse when plotted as a function 
of rescaled variable z/h(t). Therefore as turbulence de- 
velops in the domain, the level of temperature fluctua- 
tions within the mixing layer remains constant as a con- 
sequence of new fluctuations introduced by plumes enter- 
ing from unmixed regions. As Fig. U indicates, the level 
of fluctuations is larger in the viscoelastic case, as a con- 
sequence of the reduced mixing at small scales (already 
observed in Fig. [IJ . We remark that all together these 
results are consistent with the accepted phenomenology 
of viscoelastic homogeneous-isotropic turbulence where 
polymers simultaneously reduce energy at small scales 
and enhance energy contain at large scales Tol- 12|. 

The turbulent mixing layer is responsible for the huge 
enhancement of the heat exchange with respect to the 
steady conductive case. The dimensionless measure of 
the heat transport efficiency is usually given by the Nus- 
selt number Nu — (wT)hf '(kOq), the ratio between con- 
vective and conductive heat transport. For a convective 
flow in the fully developed turbulent regime, the Nusselt 
number is expected to behave as a simple scaling law with 
respect to the dimensionless temperature jump which de- 
fines the Rayleigh number Ra = Agh?/(vHi) [23J. For a 
flow in which boundary layers are irrelevant, as in our 
case, Kraichnan predicted many years ago the so-called 
ultimate state of thermal convection for which (a part 
logarithmic corrections) [23 . 24 1 



Nu = CPr^Ra 1 ' 2 
Re = DPr^^Ra 1 / 2 



(3) 



where C and D are numerical coefficients. The ultimate 
state regime has indeed recently been observed in numer- 
ical simulations of RT turbulence both in two and three 
dimensions 
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FIG. 5: Time evolution of Nusselt number Nu = (wT) /(k0 ) 
(blue, lower lines), Reynolds number Re = u rma h/p (red, 
intermediate lines) and Rayleigh number Ra = Agh 3 / (vn) 
(black, upper lines) for the Newtonian run (continuous lines) 
and the viscoelastic run (dotted lines). 
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In Fig. [5] we show the evolution of the Rayleigh number 
Ra, the Nusselt number Nu and the Reynolds number 
Re = u rms h/v as a function of time. For t > r, when 
turbulence is developed, all these dimensionless quanti- 
ties grow following dimensional predictions, i.e. Ra ~ t 6 
and Nu ~ Re ~ i 3 . Moreover, it is evident that the 
effect of polymers is to increase the values attained by 
those quantities at late time. Of course, most of this 
effect is due to the enhanced value, for the viscoelastic 
solution, of the width h(t) of the mixing layer which en- 
ters in the definition of all the quantities. As discussed 
before, another effect induced by polymers is the reduc- 
tion of small-scale turbulence in the thermal plumes, 
which leads to an additional enhancement for the heat 
flux (wT). Therefore, the Nusselt number for viscoelas- 
tic turbulence is expected to increase with respect to the 
Newtonian case when it is observed both as a function 
of time and as a function of Ra. Indeed, as shown in 
the inset of Fig. [5j both in the Newtonian and in the 
viscoelastic cases, Nu ~ Ra 1 / 2 in agreement with the 
ultimate state regime (J3j> but with different coefficients, 
C N = 0.022±0.002andCV E = 0.028±0.002 respectively, 
corresponding to an increases of 27%. 

In conclusion, we have exploited high resolution direct 
numerical simulations to investigate the effects of poly- 
mer additives on Rayleigh- Taylor turbulence. There are 
several advantages in using the present buoyancy-driven 
turbulence system. The presence of a time evolving mix- 
ing layer allows us to quantify the acceleration induced 
by polymers on a natural (nonzero) mean velocity (the 
mixing layer growth velocity) at fixed buoyancy forc- 
ing, exactly in the same spirit of usual drag reduction 
in bounded flows. The relative simple and well under- 
stood phenomenology of the heat transport (which fol- 
lows the Kraichnan's ultimate state regime) allows us to 
quantify the effects of polymers on the heat transport. 
While the former feature is specific of the present config- 
uration, the latter occurs in the bulk of the mixing region 
and therefore we conjecture that our findings hold in sit- 
uations more general than the specific setup we studied, 
as indeed a recent investigation seems to indicate [27j |. 
Moreover, RT turbulence can be realized in laboratory 
experiments and therefore our results based on numeri- 
cal simulations of primitive equations are a starting point 
of the experimental investigation of polymer additive ef- 
fects on buoyancy-driven turbulent systems. 

We thank the Cineca Supcrcomputing Center 
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sources. 
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